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ABSTRACT 



Aims. We perform a deep search for galaxies in the redshift range 6.5 < z < 7.5, to measure the evolution of the number density of 
luminous galaxies in this redshift range and derive useful constraints on the evolution of their Luminosity Function. 
Methods. We present here the first results of an ESO Large Program, that exploits the unique combination of area and sensitivity 
provided in the near-IR by the camera Hawk-I at the VLT. We have obtained two Hawk-I pointings on the GOODS South field for a 
total of ~ 32 observing hours, covering ~ 90 arcinin-. The images reach Y = 26.7 mags for the two fields. We have used public ACS 
images in the z band to select z-dropout galaxies with the colour criteria Z - Y > I, Y - J < 1.5 and Y - K < 2. The other public 
data in the UBVRIJK bands are used to reject possible low redshift interlopers. The output has been compared with extensive Monte 
Carlo simulations to quantify the observational effects of our selection criteria as well as the effects of photometric errors. 
Results. We detect 7 high quality candidates in the magnitude range K = 25.5 - 26.7. This interval samples the critical range for M, at 
z > 6 (Mi5oo - -19.5 to -21.5). After accounting for the expected incompleteness, we rule out at a 99% confidence level a luminosity 
function constant from z=6 to z=7, even including the effects of cosmic variance. For galaxies brighter than M1500 = -19.0 we derive 
a luminosity density puv = l-5;^09 ^ 10^^ erg .s"' Hz^' Mpc^^, implying a decrease by a factor 3.5 from z = 6 to z - 6.8. On the basis 
of our findings, we make predictions for the surface densities expected in future surveys, based on ULTRA- VISTA, HST-WFC3 or 
JWST-NIRCam, evaluating the best observational strategy to maximise their impact. 

Key words. Galaxies:distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: luminosity function 



1. Introduction 

The search for extremely high redshift galaxies is much more 
than an exciting exploration of the furthest frontiers of the 
Universe, although this aspect is certainly at the basis of its pop- 
ularity. Its actual astrophysical interest is tied to the constraints 
that can be set on the physical mechanisms that drove the for- 
mation and evolution of galaxies at the earliest epochs of the 
Universe. 

One important area of interest in the study of galaxies at 
z > 6 is ascertaining their role in the reionization of the Universe. 
To be fully responsible for the reionization, the density of star- 
forming galaxies at z > 7 should have been similar to that at 
z - 4-, unless there is a significant evolution in the IMF and/or 
in the dumpiness of the IGM, in the escape fraction of ion- 
ising photons or in their metallicity (see e.g. iMannucci et alj 
l2007t1Oesch et al]|2009bl and references therein). Quantifying 
their number density is therefore critical in order to con- 
strain the additional mechanisms that may be responsible for 
the re-ionizations, like Pop III dominated primordial galax- 



ies, mini -black holes or others (see e.g. IVenkatesan et al.ir2003l: 
Madau et al. 2004). 

Understanding the evolution of galaxies at high redshift is 
also very important in the broader context of galaxy evolution. 
While modelling the growth of structures of dark matter is rela- 
tively straightforward, modelling from first principles the phys- 
ical mechanisms of star formation and feedback that shaped 
galaxies across the life of the Universe is remarkably complex. 

The fundamental quantity that is cuiTently used to describe 
and quantify the galaxy population at high redshift is the UV 
Luminosity Function (LF hereafter), as derived from surveys of 
Lyman break galaxies (LBGs) at various redshifts. The evolu- 
tion of its shape and normalisation along the cosmic time pro- 
vides a clear picture of the evolution of star-forming galaxies 
in the early Universe, and an important constraint on the related 
theoretical predictions. 

Searches for LBGs have bee n extr e mely s uccessful out 
to redshift 6 (e.g. [ Steidel et a ll I1995L 119991: iBunkeretal ' 
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2004; Dickinson et al. 2004; GiavaHsco et al. 2004; Ouchi et aL 
2004; Yoshidaetal. 2006; Bouwens et al. 2003, 2006, 200^ 
Mc Lure et al.ll2009ah i.e. when the Universe was only less than 
1 Gyr old. Current samples of high-redshift galaxies (z ~ 3 - 6) 
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now contain tens of thousands of galaxies and extend to lumi- 
nosities as faint as -16 AB mag (0.01 L*). 

Despite these large samples, there is still controversy 
on how the UV LF evolves at high-redshif t. Some authors 
jSawicki & Th ompson'2006l; llwata et alfcO OT*) have argued that 
the most significant evoluti on in the UV LF happe ns at the faint 
end, others (iYoshida et all l2006t iBouwens et alj l2006i l2007h 
have foun d that the most s ignific ant evolution is at the bright 
end, while iBeckwith et alJ (l2006l) have claimed that the evolu- 
tion is similar at the bright and faint ends. Other authors have 
suggested compensating evolutions in LBG number density and 
characteristic lum inosity, resulting in a nearly constant UV lu - 
minosity density (Dickinson et al. 120041: iGiavahsco et alj|2004l) . 
The likely origins of these disc repancies are both the l arge ef- 
fect of cosmic variance (see e.g.'Trenti & Stiav ellill2008l) as well 
as systematic effects due to the different estimates of complete- 
ness level, contamination from lower redshift interlopers, vol- 
ume elements and red shift distributions in the various samples 
jStanwav et al.ll2008al) . all worsened by the known degeneracy 
among the parameters adopted to fit the LF. A clear example of 
the effect of these uncertainties is shown by the recent revision 
of the estimated slope of the UV LF at z = 2 - 3, where LBG 
samples are bigger and carefully characterised dReddv & Steidell 
I2OO9I) . 

Even within these uncertainties, it is becoming clear that the 
overall evolution of the UV LF in the redshift range z - 2 - 6 
implies a decrease in the number density of UV bright galax- 
ies (M < -20.5 ) of a factor 6-11 from z ~ 3 to z ~ 6 
(e.g. l Stanwav et ai . 2003; Shimasaku et al. 2005l: lBouwens et al] 
120061k 

On the other hand, our knowledge of the evolution be- 
yond z ~ 6 is much more scanty. Finding and study- 
ing galaxies at z - 7 to constrain the UV LF is defi- 
nitely challenging, requiring deep and wide surveys in the 
near-IR part of the spectrum. Currently, the only constraints 
come from a small sample of faint z ~ 7 - 9 candi- 
dates found in small areas within GOODS wit h deep near- 
IR J + H NICMOS and WFC3 data jB puwens & Illi ngworth 



2006"; 'Bouwens etal] 12008!: 'Oesch et al. 2009b; .Bouwens et al 



2009b; Oesch_eLalJ 0O9a; McLure et al.1 l2009bl: iBunker et al' 



2009 ). Dropout searches around lensing clusters have also 
been performed , (Richard et al. 2006, 2008; Bradley et al. 2008; 
iBouwens et al.ll200 9a; Zheng et al. 2009). The discrepant results 
of these lensing studies clearly highlight the difficulties in de- 
tecting z > 6 candidates and in re moving interlopers from such 
samples (see IBouwens et aLll2009al for a discussion). 

Further constraints come from the non-detection of bright 
candidates o f z > 6.5 galaxies in relatively wide and shallower 
observation ( Mannucci et al.ll2007l; iHenrv et al.ll2009b . 

Spectroscopic identifications of z > 6.5 LBGs are lacking 
until now, with the exception of the narrow band selected Lya 
emitter at z=6.96 by live et all (l2006l) . 

Apart from some contradictory results around lensing clus- 
ters, all the evidences suggest that the number density of UV- 
bright galaxies fades significantly at z > 6.5, amounting to a 
decrease of the volume density at the bright end o f the UV LF 
of a factor 10 -30 from z ~ 4 to z ~ 7 (Mannucci et al. 2007'; 
IBouwens et al. 2007; Stanwav et alj|200 8b; Oesch et al. 2009b) . 

It is indeed possible to use these observations to constrain 
the standard Schechter parameters of the LF, which appear to 
deviate significantly from the z = 6 ones, although the statisti- 
cal uncertainties are embarrassingly large. If one takes also into 
account the many uncertainties due to systematic error in the 



candidate selection and cosmic variance, it becomes clear that 
the evolution of the LF at z > 6 is still largely unexplored. 

To progress in this field, larger and deeper IR-based surveys 
are definitely needed. The very recent WFC3 data are provid- 
ing a dramatic advance, accessing the faint side of the LF at 
z - 7 and extending the s earches to z - 10 (Bouwe ns et al 
'2009b; Oesch et al. 20094 iMcLure et all l2009bi: iBunker et al 
2009; Yan et al. 2009). In parallel, we are conducting a search 
of z - 7 bright gala xies on wider areas using the new 
VLT IR imager Hawk-I (iPirard et al.ll2"004l: ICasah et alJl2006l; 
iKissler-Patig et al]|2008l) . which is complementary to it in many 
respects. First, we can cover areas significantly larger, albeit 
shallower than WFC3, thanks to the wide field-of view of Hawk- 

1. This will yield a statistically adequate sampling at the bright- 
est magnitudes, that is necessary to obtain an accurate estimate 
of the LF parameters. 

In addition, we will use the Y band to detect galaxies at 
z > 6.5, instead of the J band used so far Thanks to the lower 
sky background and to the extreme efficiency of Hawk-I, it is 
possible to reach the required faint limits (close to ^ 27 mags. 
AB). Such faint magnitudes can be easily reached in the J band 
only from space. More important, the shorter central wavelength 
(around \^m) of the Y filter corresponds to a lower and nar- 
rower redshift range, roughly 6.4 < z < 7.2 (see Section 4 and 
Figur^. This fair sampling is very important to verify whether 
the UV LF evolves fastly at z > 6. 

In this paper we discuss the results of the first half of our 
survey, covering a large fraction of the GOODS-S field. The pa- 
per is organised as follows: in Section 2 we present our data set 
and the multi-wavelength catalogue; in Section 3 we discuss our 
colour selection criteria and the potential interlopers affecting 
the LBG selection; in Section 4 we evaluate systematic effects 
in the light of extensive Monte Carlo simulations; in Section 5 
we present our final sample of candidate z-drop Lyman Break 
Galaxies, which is used to constrain the evolution of the z > 6 
UV luminosity function in Section 6. In Section 7 we discuss 
the implications of our findings, and of present uncertainties on 
the high-z LF, on the efficiency of future dedicated surveys. A 
summary of our methods and results is provided in Section 8. 

Throughout the whole paper, observed and rest-frame mag- 
nitudes are in the AB system, and we adopt the A-CDM concor- 
dance model (//o = 7Qkm/s/Mpc, Q.m = 0.3 and Qa = 0.7). 

2. Data 

2.1. The Data Set 

This work is based on deep y-band images obtained with Hawk- 
I, the new near-IR camera installed at the VLT. With a field of 
view of about 7.5'x7.5', pixel scale of 0.1", excellent sensitivity 
and image quality, it is currently one of the best ground-based 
instrument to search for faint, rare objects as very high redshift 
galaxies. 

We combine data collected in 2007 during the Hawk-I 
Science Verification phase with further images obtained in 2008 
through a dedicated ESO Large Programme. The images cover 
two adjacent regions of the GOODS-S field, corresponding to 
~ 85% of the deep ACS area. We name these two pointing 
GOODS 1 (southern) and GOODS2 (northern). The data ob- 
tained with the ESO LP extend the coverage of the GOODS2 
field (originally about llhrs in the Y band) to the depth of the 
GOODS 1 area. The total exposure time is 16hl5m for GOODS 1 
and 16h56m for GOODS2. The position of the two combined 
Hawk-I fields is shown in Figure [1] 
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Fig. 1. Full mosaic of the two Hawk-I images of the GOODS-S 
field. The jigsaw region shows the original GOODS ACS z-band 
image. The position of the UDF is also shown. 



The Y band images were reduced using standard techniques 
for IR data - flat fielding, sky subtraction among consecutive 
frames, and final coaddition. Particular care was taken in mask- 
ing even the faintest sources during the sky subtraction. This has 
been accomplished with a "double-pass" procedure, where the 
faintest sources have been detected in a first version of the fi- 
nally coadded images, and the whole data reduction has been 
repeated masking all objects during the sky estimate step. 

Two classes of defects are present in the Hawk-I images. One 
class consists of luminous "ripples" along the rows originating 
from bright saturated objects. These defects are due to cross-talk 
effects and have been removed by masking the whole range of 
affected rows in each image. Further problems are caused by a 
"persistence" effect, i.e. residuals left by sources which appear 
as faint objects on the next acquisition at the previous position 
(in pixel) of the bright object. These "ghosts" may appear as 
faint dropouts in the final image, unless some masking is done. 
We have decided to remove these objects from the beginning, by 
masking (in each image) the pixels where objects were detected 
in the previous one (~0. 3-1.0% of the total area in each science 
frame). This strategy turned out to be more effective and safer 
than just making a sigma-clipping during the final coaddition. 

The final images have been registered to the ACS astrometric 
solution, maintaining the original Hawk-I pixel size. 

We find fliat flie GOODS 1 and GOODS2 fields are rather 
uniform in their general properties. Through the analysis of 
bright point-like sources we determine a FWHM of 0.51 +0.01 
arcsec (- 4.8 pixels) in the GOODS 1 image and 0.49 +0.02 arc- 
sec {- 4.6 pixels) in the GOODS2 one. We computed image 
zeropoints using the standard stars observed during the same 
night of the GOODS data and at similar airmasses. Since the 
standard stars are calibrated on AOV stars in the MKO filter 
set, we have estimated the conversion to the AB system for 
the Hawk-I specific Y band filter using templates, as Ynawk-i - 



Ymko + (J MKO - Kmko)- The resulting zeropoints are Y - 26.992 
and Y = 26.998 (AB), for GOODS 1 and GOODS2 respectively. 

We have also carefully determined the r.m.s. of the coadded 
images, which is necessary to properly estimate the statistical 
meaning of detections at the faint limit. We have obtained "from 
first principles" an absolute rm.s. which fully accounts for the 
correlation in the pixel of the finally coadded image. This was 
estimated by computing the rm.s. in each individual image (us- 
ing the Poisson statistics and the instrumental gain) and propa- 
gating self-consistently this rm.s. over the whole data reduction 
process. The typical 5cr magnitude in one arcsec^ is in the range 
26.7-26.8 over more than 60% of the whole image, and > 26.2 
in 85% of the image - the rest of the images being shallower 
because of the gaps between the four Hawk-I chips. 

2.2. The Photometric Catalog - Detection 

We h ave obtained the photom etric catalog using the SExtractor 
code dBertin & Arnoutsll 19961) and the Y band as detection im- 
age. The rm.s. derived "from first principles" as described above 
is used in SExtractor as MAP_RMS and overrides the r.m.s. ob- 
tained by SExtractor from the background fluctuations. To ob- 
tain total magnitudes, we have computed both the SExtractor's 
MAG_BEST as well as aperture corrected magnitudes in 2 FWHM 
diameter (about 1"), the same apertures used for the colour esti- 
mate discussed below. We compute the aperture corrections from 
bright non-saturated stars in each field: we find corrections of 
0.364 mags and 0.322 mags in GOODS 1 and GOODS2, respec- 
tively. MAG_BEST magnitudes are more accurate for bright ob- 
jects, but become fainter than aperture-corrected ones at about 
niY - 24. As we shall show below, all z > 6.5 candidates are 
so faint that their total magnitude is estimated with aperture- 
corrected magnitudes. For resolved objects, however, aperture 
corrections based on stellar profile may underestimate the actual 
total flux. We have estimated this effect by using Lyman Break 
Galaxies with known spectro scopic redshifts 5.5 < z < 6.2 in 
the GOODS-S ACS images dVanzella et al.ll2009h . smoothing 
them to the Hawk-I PSF We find that the aperture correction is 
slightly larger, i.e. ~ 0.5. For simplicity, and given the unknown 
physical extent of our z > 6.5 galaxies, we shall adopt the nom- 
inal aperture correction based on known stars for all candidates. 
However, the effect of lost flux due to finite size of our candi- 
dates is fully accounted for by the simulations that we use to 
estimate the Luminosity Function, as we shall discuss in Section 

El 

The most critical issue is the very detection at faint limits. 
Given the extreme faintness expected for the z > 6.5 galaxies 
and their rarity, a compromise must be found between two com- 
peting goals: extending the detection at the faintest possible lev- 
els while retaining a good accuracy. In addition, the systematics 
must be understood and quantified for the final scientific analy- 
sis. We have optimised the SExtractor parameters involved in the 
detection of faint objects evaluating at the same time the possible 
contamination from spurious objects in the Y-detected catalogue 
through the analysis of a 'negative' image. To do this, we vary 
the detection parameters used by SExtractor (DETECT_MINAREA, 
DETECT_THRESH, filtering, deblending and background subtrac- 
tion parameters) in the analysis of the 'positive' image. We then 
construct a "background-subtracted negative" image and analyse 
it exactly as the positive one, with the same SExtractor parame- 
ters. The RMS image used is the same employed in the analysis 
of the 'positive' image. We finally adopt the set of parameters 
that minimises the ratio between 'negative' and 'positive' detec- 
tions at the faint end of the number counts. The final detection 
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is obtained requiring 10 contiguous pixels each at S /N > 0.727, 
corresponding to a 2.3cr detection, and restricting the analysis 
to the regions where the R.M.S is less than 1.5 times the lowest 
value. With this choice of parameters, detections on the negative 
images are negligible at F < 26.5, and less than 10% down to 
Y = 26.8 in both fields. However, a posteriori, the latter value 
overestimates the actual rate of spurious detections. Indeed, all 
spurious sources should appear as "drop-out" candidates with a 
single-band detection. On the contrary, as we will show in Sect|5] 
our two faintest candidates are both confirmed by detections in 
other IR bands. A visual inspection of the negative images shows 
that many faint 'objects' are found near bright sources: this is an 
indication that, at the faintest limits, non trivial issues concern- 
ing the subtraction of the background or a potential asymmetry 
in the noise distribution produce an overestimate of the rate of 
spurious detections. 

This procedure is also used to define the total area where 
a homogeneous catalog can be extracted. The candidates found 
in this area will be used for the evaluation of the LF. This area 
is ~ 75% of the whole image and is 89.7arcmin^. We remark 
that additional candidates satisfying our colour and S/N selection 
criteria (see Sect. [3]l are also found in the more noisy regions of 
the images, the most notable being an object falling in the UDF, 
but they will not be discussed in this paper 



2.3. The Multicolour Catalog 

We have obtained full multiwavelength photometry of all the 
objects detected on the Y band using the publicly available 
UBVRIZJHK images. ACS BVIZ images are the latest V2.0 
version released by the STSci (M. Giavalisco and the GOODS 
Team, in preparation). Each frame has been smoothed to the 
Hawk-I P SF using an appropr iate kernel obtained with Fourier 
transform (iGrazian et al.ll2006l) and registered to the Hawk-I im- 
ages. Publicly available U and R images obtained with VIMOS 
(Nonino et al. 2009) and JHK mosaics (Retzlaff et al., in prepa- 
ration) have also been registered to the Hawk-I images. These 
have not been filtered, since their PSF is always somewhat larger 
than the Y images. We then compute magnitudes in U, B, V, R, 
I, z, J, H and Ks bands running SExtractor in dual mode using 
the Y band HAWK-I image as detection image with the detec- 
tion parameters indicated above. Aperture fluxes are computed 
with the same aperture as in the Y band, and appropriate aper- 
ture corrections have been separately applied to each band. A 
comparison between ACS and ISAAC total magnitude s in our 
catalo gue and in th e GOODS-MUSIC catalogue (Grazian et all 
I2006l:lsantinietal.ll2009h shows a good agreement, apart from a 
fraction of blended objects (10% of the total) that are unresolved 
in the Y band detection but have been deblended in the ACS Z- 
detected sample of the GOODS-MUSIC catalogue. The typical 
Icr limiting magnitudes in these images, scaled to the total Y 
band flux of detected objects (i.e. estimated in the 1.2" aperture 
and corrected to total) are 29.1, 29.0, 29.1, 29.5, 28.6, 28.2, 26.6, 
26.2, 26.2 in U,B,V,R,I,Z,J,H,K respectively. 

This catalog contains self-consistent magnitudes in all 
bands. To exploit the superior image quality of the ACS images, 
we have also obtained photometry of all objects on the BVIZ 
images without smoothing them to the Hawk-I PSF, in a nar- 
rower aperture of 0.6" (an even smaller aperture would be prone 
to errors on object centering). 
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Fig. 2. Z - Y colour as a function of redshift in the Hawk-I fil- 
ter set. Shaded area shows the locus predicted by CB07 models 
with a range of metallicities, ages, dust extinction and Lyman-cc 
emission (see text for details). Red lines corresponds to extreme 
models with no Lyman-a, blue lines with a rest frame Lyman-a 
equivalent width of 300A. In the upper part, the efficiency curve 
of the two filters is shown, computed at observed wavelength of 
a Lyman-ff emission at the corresponding redshift. 



3. The Selection of z> 6.5 Galaxies 

3.1. The Colour Selection Criterion 

The selection of galaxies at z > 6.5 uses the well known "drop- 
out" or "Lyman-break" technique, with minor modifications due 
to our filter set and imaging depth. The main spectral feature 
that enables the identification of galaxies at extreme redshifts 
is the sharp drop shortward of the Lyman-a, where most of the 
photons are absorbed by the intervening HI in the intergalactic 
medium. At 6.5 < z < 7.5, this break is sampled by the large 
Z - Y colour, as shown in Figure |2] where we plot the Z - Y 
colour of galaxies at z > 6 in our filter set. At this purpose we 
have use d the models of Chariot and Bruzual 2007 (in prepara- 
tion, see lBruzuall2007a|[bl hereafter CB07) with the following 
range of free parameters: Metallicity: 0.02, 0.2 and 1 Z©; age 
from 0.01 Gyr t o the maximal age o f the Universe at a given z; 
E(B-V)=0...0.2 (ICalzetti et al.ll2000l) . Since Lyman-or emission 
has an important influence on the selection function of Lyman 
break galaxie s, as it has been proved in the case of lower redshif t 
samples (e.g. Stanwav et al.ll2008at iDow- Hvgelund et aTll2007l) . 
we explicitely take into account this effect by considering a dis- 
tribution of Ly-Qf rest-frame equivalent width in the range 0-300 
A. We have also added the intergalactic absorption using the av- 
erage evolution as in Madau ( 1995). The shaded area shows the 
region covered by the computed models. It is evident that the 
colour evolution is relatively smooth, mainly because of the ex- 
tended red tail of the Z filter which gathers a fraction of flux from 
galaxies up to z ^ 7.5. 

It is interesting to note that most of the broadening in the Z - 
Y colour distribution of the model galaxies in Figure |2]is due to 
the effect of the Lyman-c emission line, which may change the 
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Fig. 3. The distance modulus Misoo - niy as a function of red- 
shift in the F-band Hawk-I fiher. Shaded area shows the locus 
predicted by CB07 models with a range of metallicities, ages, 
dust extinction and Lyman-a emission (see text for details). 



observed colour (at a given redshift) by nearly Z - Y ^ I mags. 
The uncertainty in the intensity of this line can be translated into 
an uncertainty in the estimated redshift, which is (even ignoring 
noise effects) about Az ^ 0.2 - 0.4 in the redshift range 6.5 < 
z < 7.5. Despite this, we shall adopt a single threshold Z-Y > I 
to select candidates. 

The major difference with the standard Lyman break tech- 
nique is that the Y band that we use to detect galaxies does not 
sample the continuum around 1500 A but a region shortward of 
it, contaminated by both the IGM absorption as well as by the 
Lyman-ff emission line. Since we aim at comparing our candi- 
dates with those expected from the LF at 1500 A, this effect can 
be taken into account by computing the expected distance mod- 
ulus DM(z) - Ml 500 - my with the same template set described 
above. This is shown in Figure[3] In addition to the cosmological 
dimming with redshift (for galaxies of given M1500), the Lyman- 
a emission can produce a brightening at z > 7, rapidly counter- 
balanced at z > 7.3 by a loss of flux due to the intervening IGM 
absorption. 

All these effects will be accounted for by the Monte Carlo 
treatment that we shall discuss in the next Section. 



3.2. Interlopers 

Given the rarity of high redshift galaxies, it is mandatory to dis- 
cuss the possible contaminants in the selected colour criterion. 

Apart from z > 6 galaxies, other well known classes of ob- 
jects can display a red Z-Y colour: i) Variable objects, mostly 
due to low-intermediate redshift SN events; ii) Passively evolv- 
ing galaxies or dusty starburst galaxies at z > 1.5; iii) Galactic 
cool stars. 

We have first cross-checked each object with Z - Y > I 
against variability, by looking at images acquired at different 
epochs (1 month of delay in GOODS 1, 1 year in GOODS2). We 
identified three obvious transients in the GOODS 1 images, that 
have been removed from the following analysis. All the other 




■ 



Y-K 



I I I I I I 1 1 1 1 I 1 1 I I I I I I 1 1 I I I 1 1 1 




Fig. 4. Upper panel: Z - Y - Y-K diagram showing the ex- 
pected position of z > 6.5 galaxies (grey dots), passively evolv- 
ing galaxies and reddened starbursts (red dots) at different red- 
shifts. See text for the details of the adopted library. The red- 
shift evolution of a single representative model is also shown: 
open circles for LBGs, large filled circles for passively evolving, 
crosses for dusty starbursts. The red arrow represents the redden- 
ing vector at z = 1.5. Blue symbols mark the position of normal 
galactic stars. Lower panel: as above, for the Z-Y - Y - J colour 
plane. 



objects in our sample have a consistent photometry in the two 
epochs, at the 3cr level. 

Passively evolving galaxies or dusty starburst galaxies at 
z > 1.5 can be easily modelled with a suitable set of spec- 
tral synthesis models. We use the same CB07 library used for 
z > 5.5 galaxies to predict the colours of such objects at 1.5 < 
z < 4, using a combination of short star-formation exponential 
timescales (0.1-1 Gyrs) and ages > 1 Gyr to reproduce pas- 
sively evolving galaxies, and con stant star- forming models with 
0.5 < E(B-y) < 1.5 (adopting a ICalzetti e t al. 2000, extinction 
law) for the dusty starbursts. As shown in FiglD these galaxies 
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may have large Z-Y colours only when they also show large IR 
colour terms in the J and K bands. We note that the effect of even 
larger amount of dust extinction would be to shift the objects at 
even redder Y - J qxY - K colors. This is shown by the red ar- 
row in Fig.m which represents the reddening vector at z = 1.5. 
The slope of the reddening vector at other redshifts is very simi- 
lar. To exclude these objects, we shall adopt in the following the 
additional criteria (see Fig.|4|: 

(Z-Y) > 1.0 + {Y-K) 
(Z-Y) > 0.5 + 2.0(Y - J) 
(Y-J) < 1.5 
(Y-K) < 2.0 

It is less straightforward to exclude T-dwarfs. They are cool 
(Teff < 1500 K), low-mass stars and substellar objects that 
have been discovered and st udied in these last ten years mainly 

19991) SPS S (e.g. 



thanks to the 2MA SS (e.g. iBurgasser et al.1 

iLeggett et al.1 llOOOh. UKIDSS ( e.g. Pinfield et al.1 12008) and 



CFBDS surveys dPelorme et al.i 2008). Their infrared spectra 
are dominated by the presence of CH4 and H2O absorption 
bands and by H2 reso nant absorption (e.g. IChabrier et al.ll2005l; 
IBurgasser et alj I2OO6I) that produce a sharp break in their IR 
colours which resembles the Z-Y and Z - J colours of high 
redshift LBGs. Moreover faint dwarfs can easily be beyond the 
detection limit even in the I band because of their cold atmo- 
spheric temperatures. 

Unfortunatel y, the current uncer tainties in cool dwarfs atmo- 
spheric models dHelling et al.ll2008l) . do not allow us to compute 
the expected Y-J,Y-HotY-K colours with the reliabil- 
ity needed to distinguish T-dwarfs from LBGs. In addition, such 
distinction is unfeasible on the basis of their morphology/size in 
our ground-based Y band images, since morphological classifi- 
cation is not reliable at very low S/N ratios. 

It is however possible to estimate their expected number. 
While cool dwarfs are found also at magnitudes much brighter 
than those of z > 6 galaxies, their number density is known to 
increase towards fainter fluxes (D' Antona et al. 1999; Burgasser 
12004 ). so it is not possible to exclude that some of these objects 
contaminate our sample of high redshift candidates. The exact 
number of expected contaminants depends on the still uncer- 
tain parameters constraining the IMF and the spatial distribution 
of late type dwarfs inside the disk and the halo of the Galaxy. 
Observations have excluded that these objects can be described 
by a power-law IMF <t)(M) oc M^" as steep a s the S alpeter one 
(a = 2.35), with an upper limit at a = 1.5 (IBurg asser 2004). 
For this reason we considered a worst-case model with an IMF 
exponent a = 1.5, and a spatial distribution with height of the 
Galactic disk - 300 pc, following the T-dw arfs surface den- 
sity predictions presented by'Burgassei* ('2004') for a deep survey 
at high Galactic latitude, we obtain an estimated contamination 
of ~ 1 .7 late type dwarfs (T,,// < IISOK) for each Hawk-I point- 
ing. Larger values for the T-dwarfs scale heig ht in the disk (as the 
350 pc value proposed bv lRvan et al.ll2005l) . or an even steeper 
IMF, could slightly increase this already pessimistic estimate. 
However, we underline here that, with such a low number of ex- 
pected cool dwarfs, the main uncertainties arise from poissonian 
fluctuations in their number counts. We note that two cool dwarfs 
wi th compatible co l ours w ere indeed found in the GOODS field 
by Mannucci et alj (12007*). They bo th have z - Y > 1 in our 
catalogue: one (ID=6968 in M annucci et all (l2007h . Gl_1713 in 
our catalog having Y=24.4 5) is selected as candid ate drop-out, 
while the other (ID=4419 in lMannucci et all (l2007h . G1^130 in 



our catalog, with Y=23.3) is rejected because it has a significant 
detection also in the I band. 

Apart from these well know classes of possible interlopers 
we found that a sample of galaxies selected with the Z - Y > 1 
and Y - K < I criterion is populated also by an unexpectedly 
large number of faint contaminants showing significant detection 
in filters covering wavelengths shorter than the redshifted Lyman 
Hmit at z > 6 (U and R VIMOS; B, V, I ACS). 

We show in Fig.|5]the SED and the image in different filters 
of a typical member of this contaminating population. These ob- 
jects show a dip in the observed flux in the Z and often in the 
/ band, with a rising continuum in the bluer bands and a large 
Z-Y colour. Such a spectral energy distribution cannot be re- 
produced by a straightforward application of the CB07 models. 
The analysis of the nature of these objects is beyond the scope 
of the present paper and will be caiTied out separately. We only 
note that all of them are relatively faint (mostly my > 25). A vi- 
sual inspection of their ACS images revealed that these objects 
show a variety of morphologies, at least those detected at good 
S/N. In many cases they are clearly extended, thus indicating 
their extragalactic nature, although other cases present a point- 
like morphology. One possibility, especially for those undetected 
in JHK, is that they are faint galaxies with a very blue contin- 
uum whose SED is altered by the presence of strong emission 
lines such as in un-obscured AGN, o r in star f orming galaxies 
like the blue compact dwarf galaxies dlzotov et a l. 2 004, 2007!) 
or the Ultra Strong Emission Line Galaxies (USELs, iHu et al.l 
200^. 

Although standard prescriptions for the selection of high 
redshift LBGs already exclude sources detected at short wave- 
lengths, given the rarity of z > 6 candidates we took particu- 
lar attention in tailoring reliable criterion to separate our sam- 
ple of high redshift candidates from these unexpected lower red- 
shift objects. After looking at the colour distributions of these 
contaminants and at the distribution of the detected S/N in the 
UBVRI bands, we decided to adopt a selection criterion even 
more conservative that the one typically adopted by previous sur- 
veys. We require the S/N from high-z candidates to be < Ios/n 
in all UBVRI bands and < Ios/n in at least four of them, os /n is 
the estimated rm.s. of the S/N distribution in each band, as de- 
scribed in the next Section. In the case of the ACS BVI images 
we used S/N ratios, and relevant cts/n, measured in a smaller 
aperture (0.6", see Sect. 12.31 ) to better exploit their higher res- 
olution. To check whether our selection criteria are effective in 
excluding this population of faint contaminants from the LBG 
sample, we have also analysed the d eep (~ 30 mags AB at Icr) 
public images of the UDF area (Be ckwithet"an 120061) . which 
is also covered by our shallower data-set (see Fig. [T]). Those 
objects that, on the basis of the conservative criteria discussed 
above, are individuated as contaminants in our catalogue, are all 
confirmed, at an higher significancy, to have emission in one or 
more of the blue bands. On the other hand, we effectively select 
as an hig h redshift candid ate object UDF-387-1 125 already dis- 
cussed in iBouwens et al.l (l2004j) and rece ntly confirmed in th e 
analysis of the very de e p WFC3 images byfeunker et aU ( 2009 1; 
iMcLure et all (l2009bl) : lOesch et all (l2009al) : lYan et alJ (l2009l) . 
Thus we can state that deep optical images as the ones analysed 
here (~ 29.0 mags AB at Icr) are sufficient in separating this 
class of contaminants from more reliable high redshift candi- 
dates, once conservative selection criteria are adopted. However, 
we caution that an accurate, dedicated, spectroscopic analysis of 
both LBGs and contaminants will be necessary to determine the 
real impact of these objects in LBG searches, since their physical 
nature is still undetermined. 
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Fig. 5. SED and thumbnail for one of the interlopers with Z- Y > 
1 and detection at shorter wavelengths. 



4. Simulating tlie Systematic Effects 

Although the colour criteria are formally very clear-cut, they are 
in practice applied on very faint objects, typically close to the 
limiting depth of the images. At these limits, systematics may 
significantly affect their detection and the accurate estimate of 
their large colour terms. To fully evaluate the involved uncer- 
tainties we have performed extensive imaging simulations, by 
which we have quantified the systematic effects in the object de- 
tection, in the measure of their total magnitude as well as in the 
measure of large colour terms in our images. The output of these 
simulations will be used in the proper estimate of the LF. 

At this purpose, we first use the synthetic libraries described 
above to produce a large set of simulated galaxies with expected 
magnitudes in the UBVRIZY filter set in the redshift range 
5.5 < z < 8. Objects are normalised in the range F = 24 - 27.5, 
following the expected magnitude distribution arising from a LF 
with index a - -1.71, These galaxies are placed at random po- 
sitions of the GOODS 1 and GOODS2 F-band images, and cata- 
logs are extracted exactly as in the original frames. 

As expected, the output of this exercise depends critically 
on the assumed morphology. In our case, the availability of 
deep z-band ACS images of confirmed z - 5.5-6 LBGs pro- 
vides the most natural templates, avoiding further assumptions. 
Neglecting possible size evolution from z = 6 to z = 7, we have 
used the four brightest z = 5.5 - 6 LBGs observed with ACS in 
GOODS, both convolved with the GOODS 1/2 PSFs. 

To avoid an excessive and unphysical crowding in the sim- 
ulated images, we have included only 200 objects of the same 
flux and morphology each time, after masking the regions of the 
images where real objects have been detected (this corresponds 
to 6% of the total area). We repeated the simulation until a total 



I 
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Fig. 6. Upper pane tFTaction of detected objects as a function 
of their input magnitude, as estimated from the simulations 
described in the text. Lower panel: From the same simula- 
tions, contour levels of the conditional probability P{Ym, (Z - 
Y)m\Y, (Z-Y)) that a galaxy with given Y band flux and Z-Y 
colour is detected with a measured Ym flux and with a measured 
(Z - Y)m colour. Only four cases are shown, with input values 
marked by large crosses. Contours are drawn at 0.05, 0.15, 0.3, 
0.5 and 0.8 times the peak value of each distribution. 



of 10^ objects were tested for each field and morphology, for a 
total of 8 X lO*" artificial objects. 

These simulations can be used to estimate the impact of dif- 
ferent systematics in the various steps of the analysis. First, they 
can be used to estimate the completeness in the detection pro- 
cedure. This is shown in the upper panel of Figure |6] where we 
plot the fraction of detected versus input objects as a function 
of the input magnitude. The detection is close to 100% down to 
Y ^ 26.2, and fades to 30% al Y ^ 26.7 (where we find our 
fainter candidates. We remind that an additional 6% in the in- 
completeness is due to the masking done to avoid bright objects. 

In addition, they can be used to evaluate the uncertainties in 
the estimate of the colour criteria that we shall exploit to detect 
z > 6.5 candidates. In particular, the Z-Y colour is critical fea- 
ture to identify z > 6.5 objects. To visualise the eff'ect of noise, 
we plot in the lower panel of Figure|6]the conditional probability 
P{Ym, {Z-Y)m I Y, (Z-Y) ) that a galaxy with given Y band flux 
and Z-Y colour is detected with a measured Ym flux and with a 
measured {Z-Y)m colour. Obviously, a fraction of galaxies will 
not be detected at all in the Y band, as discussed above. 

An inspection of Figure |6] shows the basic feature of the sys- 
tematics acting on our images. At relatively high S/N (Y ^ 25, 
Z ^ 25), the recovered magnitudes and colours have a narrow 
scatter, cry - cr, ^ 0.07 mags which increases at lower fluxes 
( o-y ^ ^ 0.2 mags at F ^ 26, Z 26). At larger Z-Y 
colours asymmetries becomes evident as the Z magnitude ap- 
proaches the very detection limit. At the faintest limit, {Y ^ 26.8 
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Fig. 7. Upper panel: Fraction of F-detected objects passing our 
criteria {Z - Y > 1 ; non-detection blueward of Z) as a function 
of redshift, for two limiting magnitudes, Y - 26.3 (solid line) 
and Y - 26.7 ( dashed Une). I nput objects are extracted assuming 
a constant LF jMcLure et al., ,2009a) from z - 6. Lower panel: 
Redshift distribution (per unit redshift) for the same sample. 



and Z -Y ^ 1), about 30% of the objects detected in Y become 
too faint to be detected in the Z band image. 

The simulations in Figure |6] also show that a 0.15 mags off- 
set exists on average between the input and the recovered mag- 
nitude. As discussed above, this is due to the adopted aperture 
correction, which is computed on unresolved stellar profiles, in- 
stead of the LBG profile adopted in the simulations. 

Finally, this set of simulations is used to estimate the system- 
atic effects acting when we use colours at shorter wavelengths, 
i.e. in the UBVRI bands. In these cases, because of the large 
IGM and internal HI absorption, the expected flux in these bands 
for z > 6.5 galaxies is far below the detection threshold, or 
even null. For this reason a stringent limit on the measured flux 
in these bands is adopted to remove lower redshift interlopers. 
However, the S/N estimated by SExtractor may be a poor repre- 
sentation of the actual photometric scatter at low fluxes, due to 
a combination of factors, such as uncertainties in the estimate of 
the local background, underestimates of the true r.m.s., or chance 
superposition of faint blue galaxies along the line of sight. 

To account for these effects, we have measured the resulting 
signal-to-noise SN in the UBVRI images for each simulated ob- 
jects inserted in the Y one, which should be zero on average. It 
turns out that the actual distribution of the S/N ratios is wider 
than the one obtained with SExtractor, which is computed scal- 
ing the input weight image. With this set of simulations we thus 
estimate the "effective" r.m.s. os/n, i e. the rm.s. of the signal- 
to-noise distribution in each of the 5 images, which is typically 
about 1.5. Even taking into account this wider distribution, we 
also find that the tails of the S/N distribution (S/N > Ios/n) 
contain more objects than in the case of a pure Gaussian distri- 
bution. 

As mentioned above, we will use the estimated os/n in 
all UBVRI bands, requesting that high-z candidates have flux 
< IcTsiN in all UBVRI bands and < Ictsin in at least four 
of them. With our simulation, we estimate that the fraction of 



true high-z galaxies lost because of this strict criterion is about 
~ 30%. This effect will also be taken into account in our estimate 
of the LF 

The output of this exercise can be summarised in Figure |7] 
The final effect of noise on the Z - F > 1 colour is shown in 
the upper panel, where we plot the fraction of simulated objects 
detected in Y and with measured Z - Y > 1 as a function of 
redshift, for two limiting magnitudes, Y - 26.3 (solid line) and 
Y - 26.7 (dashed line). We remind that, since the objects are 
extracted from a steep LF, most of them are close to the limiting 
magnitude. The Z-F > 1 colour and the requested non-detection 
blueward of Z provide an effective cut below z - 6.3-6.4, which 
is nearly total at z < 6, even including the effect of noise. At 
z > 7, the fraction of detected objects is less than 100% because 
of the combined effects of incompleteness in the detection and of 
the requested non-detection blueward of Z, which rejects some 
genuine high-z galaxies because of photometric scatter. 

The most important output is shown in the lower panel of 
Figure|2l where we plot the expected redshift distribution of our 
sample, as extracted from the simulations. In this case, counts are 
normalised assuming that the LF remains c o nstant from z = 6 
with the values provided by iMcLure et al.l (l2009al) . Since we 
populate the input catalog following a LF with a realistic slope, 
we expect that this redshift distribution is a good estimate of 
the redshift selection function of our survey. The low redshift 
cut-off is due to the Z - F > 1 colour and the requested non- 
detection blueward of Z, as discussed above. At higher redshift, 
the cutoff is due to the decreasing number of expected galaxies 
in a magnitude-selected sample - as shown in Figure [3l objects 
of given L become rapidly fainter at z > 7.3. In this case, the 
high redshift tails of the distribution is sensitive to the evolution 
of the LF: in case of a mar ked drop of luminou s objects at high 
redshift, with respect to the McLure et al.l (12009 a ) LF, the distri- 
bution will be skewed toward the lower redshift boundary. This 
plot clearly shows the basic feature of our F-band selected ap- 
proach, which allows to sample the LBG LF in a well defined 
narrow redshift range. We note that this redshift range is nar- 
rower than the corresponding one in z - / selected samples . 

5. Detected z> 6.5 Galaxies 

Based on the results of the previous sections, we finally select 
the z > 6.5 candidates adopting the following criteria: 

1) We include only objects detected in the deepest regions of the 
images (r.m.s. not larger than 1 .5x its minimum value) with total 
magnitude Y< 26.7; 

2) We require: 

(Z-F) > 1.0 

(Z-F) > 1.0 + {Y-K) 

(Z-Y) > 0.5 + 2.0(F - J) 

(Y-J) < 1.5 

(Y-K) < 2.0 



3) We require that the S/N (as measured in units of the 
simulation-calibrated crs /n) is less than 2 in all UBVRI bands, 
and above 1 in one of these bands only. 

Figure [8] shows the position of all objects of our sample in 
the F - 7 vs Z - F plane, with the final candidates and the two 
known brown dwarfs shown highlighted. 

In practice, the limit on the Y-J < 1.5 and Y-K < 2 colours 
are ineffective. Indeed, all the objects with Z - F > 1 are very 
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Fig. 9. Thumbnails showing the images of the 7 selected high redshift candidates in the different observed bands. 



blue in the near-IR, having Y - J < 1.2 and Y - K < 1.7, and 
90% of them have F - 7 < 0.5 and F - < 1. 

As a result, the critical requirement to select "bona-fide" 
high-z candidates is the non-detection in all the bands blueward 
of the Z filter. In our catalog we detect 1 39 objects with Z-Y > I 
and Y < 26.8, only 9 of which pass the additional criteria on the 
non-detection in the UBVRI bands. 

Most of our candidates are marginally detected at very low 
S /N - 2 in the Z and/or J images, which makes mandatory to 
critically assess their reality. The first and most obvious check 
to be done is with the available NICMOS FllOW and F160W 
images in the GOODS area. Unfortunately, these images do not 
cover with the required depth the whole Hawk-I pointings. Out 



of 4 of our candidates falling on suitably deep images, we clearly 
detect 3 of them with Fl lOW and F160W magnitudes consistent 
with our Y one. The fourth one is undetected in the F160W with 
mpibow - my > 1 (FllOW is not available). As a consequence 
we remove this object from our sample. 

A similar check can be done on the IRAC images. Taking 
into account the actual depth of the IRAC images, a detection 
is not definitely required to confirm the existence of these candi- 
dates. Indeed, using the same synthetic libraries described above, 
we expect colours my -m/sAc in the range -2/ +2 mags, primar- 
ily depending on the age, mass and dust content of the galaxies. 
Since the detection limit in IRAC does not exceed 26.5, grossly 
the same AB limit of our candidates, we can expect that only 
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Fig. 8. Position of the high-z candidates in the Z - F vs F - 7 
colour plane (black large dots). Upper limits are computed at 
the Icr level. The shaded area shows the location of the expected 
colours of CB07 models, at z > 6.5. The faint gray dots show the 
rest of the sample. Objects fulfilling the selection criterion are 
rejected because of their detection in at least one of the UBVRI 
images (see text for details). The two Brown Dwarfs akeady 
identified by Mannucci et al. are shown as stars. The one in the 
upper right is selected as a high-z candidate on the basis of its 
UBVRIZY photometry, but removed from the final sample. 



a fraction of these galaxies are detected in the IRAC images. 
Indeed, 3 (out of 8) objects are clearly detected in the IRAC 3.6jj 
channel (IRAC36 hereafter): Gl_1921, G2_1713 and the object 
G 1^631 already detected in NICMOS. Further three objects 
(G2_1408, G2_2370 and G2_6173) are mai-ginally detected due 
to blending with nearby sources. We note that two of the IRAC 
detected objects are the faintest ones in our sample (Gl_1921 
and Gl_2631), so the IRAC detection confirms their reality. 

We computed aperture photometry (corrected to total) on 
both the 3.6/i and the 4.5yL( IRAC image. None of our IRAC de- 
tected objetcs shows a (3.6 - 4.5) index larger than ~ 0.4, in 
agreement with the colours expected for high redshift LBGs. On 
the other hand, we note that previous studies have found that cool 
dwarf stars show a broad range in the 4.5/i luminosity, yield- 
ing values in the 3.6 - 4.5 colour index that can be muc h larger 
than those of our objects ( up to 2.0 magnitudes, see IPatten et al.l 
120061: iHelling et al.ll2008l) . although lower colour indexes have 
also been observed for these objects. 

As already described in Sect. 3.2, one of our candi- 
dat es was already detected as a potential high-z candidate 
by iMannucci et al.l (l2007l) and associated to a Galactic brown 
dwarf. 

We show in Figure |8] the p osition of this and of the other 
Galactic brown dwarf found by IMannucci et al] (l2007h . as well 
as of the other 7 candidates. It is interesting to note that, given the 
high S/N that we obtain on these relatively bright brown dwarfs, 
the measured Y - J colour is larger than that expected for a 
Lyman-break gala xy at z > 6.5. This an alysis further supports 
the conclusions of iMannucci et al.l (l2007b . therefore we remove 
Gl_1713 from our sample of "bona-fide" high-z candidates. 



In the attempt to identify further brown dwarfs, we have also 
investigated the morphology of the candidates in the NICMOS 
images, where available. Object G2_2370 has a full width of 
about 0.45", close to the 0.4 value in NICMOS, but with a 
non-null ellipticity 0.3. Object G2_1408 is clearly resolved, with 
FWHM^ 0.65. Considering also that several LBGs with spec- 
troscopic redshifts in GOODS have a NICMOS PSF consistent 
with the stellar one, we shall keep both of them in our sample. 

After removing from our sample the IMannucci et all (|2007|) 
brown dwarf and the object undetected in deep NICMOS point- 
ings, we are left with 7 candidates: their position and properties 
are given in Table 1 . The relevant thumbnails are given in Figure 

121 

We performed a stacking of all the thumbnails in the 
UBVRIZYJHK images. We confirm that objects are undetected 
in the UBVRI images, are detected in the J band average image 
and exhibit an average colour Z - F ^ 2. The resulting SED 
provides an excellent photometric redshift at z = 6.8. Relevant 
thumbnails and SED are shown in Figure [TOl 

We checked if any of the previously published, NICMOS de- 
tected, z > 6 candidate s in the UDF area (B ouwens et al. 200^ 
iBouwens & IlUngworthl 120061: lOesch et all l2009bl) are present 
within our sample. We fi nd that our can d idate G2_1408 is the 
object UDF-387-1125 in iBouwens et alJ (l2004 . All the other 
objects are not detected with our extraction parameters. We no- 
tice however that most of these objects have faint magnitudes in 
the JllO and H160 NICMOS bands (~ 26.8 - 27.0 AB), while 
we adopted conservative selection criteria in order to have a low 
number of spurious detections at the faintest magnitudes in our 
sample (Y - 26.8, see Sect. 12. 2t , so it is not surprising that we 
are missing objects detected in the very deep NICMOS point- 
ings over the UDF. We note that the colours of all the selected 
objects are consistent with small values of F- 7 (or Y-K) plane, 
as shown in Figure [8] Since larger value of F - 7 (even within 
our threshold Y - J < 1.5) are typical of galaxies at z > 7.2, 
this confirms our expectations that the redshift distribution is in 
practice limited at this value. 

We also checked other recent results on the GOODS area 



pre sented bvlHickev et al. (2009) and by Wilkins et al. (2009). 

iHickev et alj (120091) exploit part of the same Hawk-I data-set 
presented here but reach a shallower magnitude limit (Y~25.7- 
25.9). Their candidate ID 9697 is our candidate G2_2370, while 
their object ID 9136 is rejected in our sample beacuse of S/N 
in the I band slightly above our threshold. Their other two can- 
didates show significant detection in the bluer bands and so are 
probably lower redshift interlopers as discussed also by the au- 
thorS; 

IWilkins et all ([2009) present an analysis of the WFC3 ERS 
data covering the northern portion of our "GOODS2" field. Their 
candidate #4 is our object G2_6173, although they measure a 
fainter Y magnitude. Of their two candidates having Y<26.8 we 
reject object #1 beacuse of non negligible S/N in the I band, 
while object #2 is too close to a bright (Y~20) source to be 
effectively de-blended in our ground-based images and it is not 
present in our catalogue. 

We remark that we estimated that a 30% of real sources are 
missed with the strict rejection criteria in the blue bands adopted 
in the present analysis . This is consistent with object ID 9136 
in I Hickev et al.l (|2009|) and object #1 in IWilkins et all (|2009|) 
being real high redshift galaxies whose high I-band S/N in our 
catalogue is due to photometric scatter, although their non negli- 
gible S/N in the I band might be due to a redshift near the lower 
limit of the redshift selection window. We remark that potential 
missed detections beacuse of photometric scatter, or beacuse of 
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Table 1. Hawk-I z-drop candidates in Goods-S 



ID 


R.A. (deg) 


DEC. (deg) 


Y 


Z-Y 


Gl_149 


53.157038 


-27.930707 


26.06 


1.60 


G1.1921 


53.19538 


-27.835421 


26.71 


1.02 


G 1.2631 


53.106022 


-27.848123 


26.64 


1.06 


G2_1408 


53.177382 


-27.782416 


26.37 


>2.1 


G2_2370 


53.094421 


-27.716847 


25.56 


1.68 


G2_4034 


53.150019 


-27.744914 


26.35 


>2.1 
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Fig. 11. Y and Z-Y distributions of the high redshift candidates 
and of our best-fit evolving LF at z > 6 (black solid line), com- 
pared with the expected va l ues fro m recent estimates of the LF 
at z ~ 6 by Bouwens et al. (2007) (B07, green dot-dashed line) 
and|Mc Lure et aL^ (2009a) (M09, red dashed line), as computed 
according to the simulations described in Sect. |6] 



Fig. 10. L'p/^er.- Thumbnails showing the stacked images of the 
7 selected high redshift candidates in the observed bands shown 
in the legends. Lower, resulting SED, with relevant photometric 
redshift at z = 6.8 



objects falling near other bright sources, are fully taken into ac- 
count in the LF estimate presented in the next section based on 
extensive imaging simulations. 

We finally checked whether the re is any X-ray e mission de- 
tectable in the deep X ray images jLuo et al.ll2008l) . No source 
was individually detected. In the stacked image, the 90% count 
rate (0.5-2keV) limit is 2.41 x 10"^ counts/s, for a total exposure 
time of 1.2xl0^s. The resulting Luminosity is < 7.8x10"*^ ergs/s 
in the band from 2 to 10 keV (assuming a power-law spectrum 
with a energy index -0.4). Assuming a bolometric correction fac- 
tor of 20 and a Eddington-limited accretion rate, we derive a 
limit on the black hole mass of about < IO^Mq. 



6. The Evolution of the LF 

To estimate from our detections the most likely LF we use an 
approach that fully accounts for the exp ected systematics i n 
the detection process, a s done by, e. g., Bouw ens et al.l (l2007l) : 
iMannucci et all (|2007|) : iMcLure et ai., (2009a) . We refer in the 
following to the luminosity L as measured at 1500A rest-frame. 

First, we assume that the LF can be described by the usual 
Schec hter function with parameters 0, a and M, (Schechter 
Il976l) . Unfortunately, we are unable to constrain the slope of 
the LF, since our faint limit is close to the expected value of 
the characteristic luminosity M». For this reason, our results are 
grossly insensitive to the value of the slope index a, that we shall 
fix to th e value a - -1.71 of the z ~ 6 LF by McLure et al.l 
( l2009al) . close to a = -1.74, used bv lBouwens et all ( l2007h . We 
have explicitly tested this assumption by fixing a to significantly 



different values (a — -1.4, -2.0), without finding significant dif- 
ferences. 

We will assume that the other two parameters evolve in red- 
shift from their value at lower z - Zo value: 

log{(l)(z)) = logiipizo)) + dlog(cp)ldz ■ {z - Zo) 

M.(z) = M.(zo) + m: ■ (z - zo) 

We not e that our para metrisation is analogous to the one 
adopted bv lBouwens et al.1 (|2008). While iBou wens et al.l(l2008h 
used this parametrisation to fit their values at different redshifts 
(in the range 4 < z < 10), we explicitly vary the Schechter pa- 
rameters within our redshift range 6.4 < z < 7.2, which might 
be important in case of a strong evolution. We have also tested a 
purely linear evolution of (p, which turned out to produce consis- 
tent results. 

In principle, all four parameters could be left free in the min- 
imisation process. However, given the size and depth of our sam- 
ple, we shall assume for M,(zo) and 0(zo) the observed values at 
slighter lower redshifts z - 6 and evaluate the evolutionary terms 
M'^ and dlog((p)/dz alone. 

During the minimisation process, for any given value of the 
free parameters, we Monte Carlo extract a sample of objects 
with redshift z, 1500A rest-frame absolute magnitude M, ran- 
domly adding diff'erent E(B-V) and metallicities as obtained by 
the CB07 models shown in Figure |3] and Figure |2] We include 
Lya emission with a Gaussian distribution with standard devia- 
tion 30A. All these galaxies are extracted from the (larger) simu- 
lations described in Sect. |4] This way, we can convert their rest- 
frame Muv into observer-frame magnitudes in the same bands 
used to in our observed sample, taking into account all the un- 
certainties involved in the observations discussed in Sect. H) de- 
tection completeness, photometric scatter and random fluctua- 
tions in the S/N measure due to overlapping interlopers or other 
effects. 
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Fig. 12. contour levels for the dlog((p)/dz. Ml parameters de- 
rived for the Schechter-like LF. The lower and left axis refer to 
the evolutionary terms M'^ and dlog{(p)/dz (see text for details), 
the upper and right axis refer to the M, and <p values at the me- 
dian redshift estimated for our sample (z=6.8). 
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Fig. 13. 68% (dotted line) and 95% (continuous line);^'^ contour 
levels for the value of M, and (p atz - 6.8 as estimated from our 
sample of Y-detected LBGs. Blue, green, red and magenta con- 
tours are the 68% and 95% likelihood intervals on the Schechter 
parameters as estimated by Bo uwens et al.. (.2007.) at redshifts 
3.8, 5.0, 5.9 and 7.4 respectively. 



The distributions of magnitudes and colours for each Monte 
Carlo simulation are scaled to the observed area in the GOODS- 
S field (after excluding the fraction of area lost because of the 
presence of lower redshift objects) and compared to the observed 
ones with a maximum likelihood test under the assumption of 
simple poissonian statistics. 

To provide a visual example of the significance of the re- 
sults, we show in Figure[TT]the expected number counts for three 
cases of the evolution of the LF, comparing them with our ob- 
served counts. It is immediately appreciated that, barring for the 
moment large fluctuations due to cosmic variance, the number 
counts of our candidates imply a strong evolution of the LF in 
the relatively short cosmic time elapsed from z ^ 6 to z = 6.8. 

To constrain the evolution of the LF, for each simulated pop- 
ulation, and for each of the two distributions, we build the like- 
lihood function X. 



(1) 



where Nohsj is the observed number of sources in the magnitude 
(colour) interval /, Ne^pj is the expected number of sources in the 
same magnitude (colour) interval, and 11, is the product symbol. 
The final likelihood assigned to each given model is the prod- 
uct of the likelihoods computed for the magnitude and colour 
distributions separately. We iterate this procedure for a grid of 
(dlog{(p)/dz, M'J values. 

The colour plot in Fig.[T2]shows the 68%, 95% and 99% like- 
lihood intervals on the evolutionary terms M'^ and dlog{(p)ldz 
(left and bottom axes) and for the resulting Schechter parame- 
ters at the median redshift z=6.8 of our sample (top and right 
axes). In the same plot, the colour code refers to the;t'^ distribu- 
tion obtained under the usual assumption - -2.0 • ln{£). It 
is evident that the absence of any evolution in both parameters 
(dlog((p)/dz - and M'^ - dM^/dz - 0) with respect to the best- 
fit values at z = 6 is ruled out at > 99% confidence level. This is 
shown more clearly by Fig. [T3] where our 68% and 95% likeli- 
hood intervals at z=6.8 are displayed together with the likelihood 
contou rs computed in a self-consistent way by Bouwens et al] 
(120071) at redshifts 3.8, 5.0, 5.9 and 7.4. 

The formal maximum of the maximum likelihood lays at 
dlog(<p)/dz = -0.89 + 0.21 and M' = dMJdz = -0.25 + 0.56 
(errors include the effect of cosmic variance, see below). At face 
value, our best fit model indicates mostly a decrease in the nor- 
malisation factor <p when compared to the best fit parameters at 
Z ~ 6 (Fig.[T3]l. However, the relatively small number of galax- 
ies that we have in our fields does not allow us to put tight con- 
straints on the exact combination of Schechter parameters: as 
shown by the elongated shape of the contour levels in Fig.[T3]the 
two parameters are highly degenerate. We indeed find that the 
position of the maximum is poorly constrained, and somewhat 
dependent on the details of the simulations, like step or number 
of simulated galaxies. We therefore provide its value only for 
reference, and will focus on the overall range of allowed values 
at the Icr contour or on integrated quantites. 

The effects of cosmic variance are a significant concern in 
our case, since all our data come from a single pointing which 
includes most of the previous NICMOS-based surveys. Indeed, 
cosmic variance ef fects are not independ ent from those poten- 
tially affecting the iBouwens et al.l (l2008l) results. The cleanest 
way to get rid of these effects is the analysis of several inde- 
pendent fields applying the same selection criteria. For this rea- 
son, our survey will continue by covering two completely dis- 
joint fields with the same setup. For the moment, we can evaluate 
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Table 2. Stepwise determination of the UV LF at z~6.8 





<p {\Q-*Mpc-' mag-') 


-22.0 


02+"-^^ 
-Q-Q3 


-21.0 


19+"'*' 


-20.0 


2 20+'-'" 



the possible impact of cosmic variance by measuring the relative 
variance within 200 sampl es bootstrapped from t he Mi llennium 
Simulation presented by iKitzbichler & White! (l2007l) . that is 
considered to accurately reprodu ce clusterin g properties even at 
very high redshift (e.g. Overzier et al.ll2009l) . We use an area as 
large as our Hawk-I data set and we apply a corresponding pho- 
tometric selection criteria on galaxies at 6.4 < z < 7.2, with- 
out any constraint on the distribution of host halos. We estimate 
that a cosmic variance of ~ 30% affects the number counts of 
z-drop LBGs in our ~ 90arcmin^ area. We have estimated the 
resulting effect on the (p and M» parameters changing by +30% 
the observed number counts, and finding the relevant best-fit pa- 
rameters. These offsets were added in quadrature to the Hessian 
errors computed around the best fit value. 

In this way we obtain a best fit interval for the evolution- 
ary terms along with their associated uncertainty which in- 
cludes both Poisson and cosmic variance noise: dlog{(p)/dz - 
-0.59 ± 0.25 and M' = dMJdz = 0.13 ± 0.51, resulting in 
4>* = 0.35+° J^10-3Mpc-3 and M*(1500A) = -20.24 + 0.45 at 
z=6.8. 

Fortunately, the degeneracy in the M, and (f> best-fit values is 
not reflected in a comparable uncertainty in the number density 
of bright galaxies, i.e. in the integral of the LF up to M». Indeed, 
because of the correlation between the acceptable values of M» 
and 0, an increase in the normalisation {<p) is compensated by a 
decrease in the characteristic luminosity, such that the integral 
up to M ^ M, remains grossly constant. Integrating the best-fit 
UV LF up to M = -19.0 we obtain an UV luminosity density 
puv - l-5^Q9l0^^erg s"' Hz^^ Mpc^^. The errors are computed 
integrating all the UV LF that are acceptable at the 68% level. 

For comparison, the integral of the z=6 UV LF of 
iMcLure et al.l (l2009ah up to the same magnitude limit yields a 
Puv = 5. 6'^2^1Q^^ erg i"' . Our estimate thus implies a drop of a 
factor ~ 3.5 in the UV luminosity density from z=6 to z=6.8. 

We can convert this value in a star formation rate density 
following the standard formula by Madau et al. ( 1998) and ap- 
plying the extinction correction of iMeurer et alj ([1999) (con- 
sidering an average UV slope /3 - -2.0). We obtain SFRD - 
2.r}^^\Q-^Msunyr-' Mpc-\ 

To provide a more straightforward comparison to other re- 
sults in the literature, we have also computed the luminosity 
function for our z-dr opout sample in a stepwise form (see, e.g. 
iBouwens et al.ll2008l for a discussion of this method). Briefly, 
the stepwise method assumes that the rest-frame luminosity 
function of galaxies can be approximated by a binned distribu- 
tion, where the number density 0, in each bin is a free parameter: 
this non-parametric approach allows us to constrain the number 
density of galaxies at different magnitudes without assuming a 
Schechter-like shape. 

We have assumed that the LF is made of three bins in the 
interval -22.5 < M1500 < -19.5, corresponding to the range 
sampled by our observations. We also assume that galaxies are 
uniformly distributed within the bins, with number densities 0, 
to be determined. We exploit the same set of simulations de- 
scribed in Sect. |4] to compute the distribution of observed mag- 
nitudes originated from each bin, scaled to the observed area in 




-22 -20 -18 

M1500 



Fig. 14. Number densities in three rest-frame magnitude in- 
tervals estimated for our Hawk-I data-set in a stepwise form 
(black circles and error bars), along with the results by 
Bouwens et al. (2008) for NICMOS detected objects (red fifled 
squ ares), Ouc hi et al. (2009) (SUBARU, blue empty squares) 
and lOesch et a l. (2009a) (WFC3, magenta empty circles). The 
black solid line is our best-fit LF discussed in the text. For a 
compariso n we s how the recent determinations of the LF at 
z ~ 6 by iBouwen s et al. (2007) (B07, green dot-dashed fine) 
and lMcLure et all (12009a) (M09, blue dashed fine). 

the GOODS-S field. We then find the combination of 0, fliat best 
reproduces the magnitude distribution of our observed objects 
with a simple minimization. 

The results are reported in Table |2] and are displayed in 
Fig. [T4]along with recent results from the literature and our best- 
fit Schechter LF discussed above. 

Our two independent determinations of the LF at z~ 6.8 
(stepwise and maximum likelihood) are in perfect agreement. 
Th e agreement at the br ight end with the den sities estimated 
bv lBouwens et al.1 (l2008h and lOuchi etal] (12009 ') is remarkable, 
considering the widely different data-sets and selection tech- 
niques used and the large uncertainties given by the low number 
of objects in our sample. The disagree ment at the fa i nt end b e- 
tween our best-fit LF and the points by lOesch et"an (l2009ah is 
probably the result of the fixed slope a - -1.71 we had to as- 
sume for our maximum likelihood test, as discussed above, and 
it is consistent with a steepening of the faint end, as suggested 
by the authors. 

7. Predictions on future surveys 

It is interesting to use the values of the Schechter parameters 
found in the previous section to estimate the number of detec- 
tions that may be obtained with future surveys. Clearly, the large 
errors in our determination make this exercise uncertain, but nev- 
ertheless useful to design future surveys from ground and space. 

From the ground, the relative efficiency of the Y band al- 
lows us to cover significantly large areas of the sky. This is use- 
ful both to beat cosmic variance as well as to identify potential 
spectroscopic targets. Using our best fit values of dlog{(p)/dz, i» 
(and considering the two values a'=-1.7 and q'=-1.9) in Figure[T5] 
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Fig. 15. Expected number of 6.4 < z < 7.2 galaxies as a 
function of area and limiting magnitude for present and future 
ground-based surveys and for deep pointings with HST-WFC3 
and JWST-NIRCam (see text for details). Black curves show the 
combination of area and limiting magnitudes necessary to col- 
lect 1, 10, 100 or 1000 galaxies in this redshift range. Estimates 
are based on our best- fit values of dlog{(p)ldz, M[ and consider- 
ing the two values a--\J and a--\.9 (upper and lower curves, 
respectively). Note that incompleteness and noise effects are not 
included in the computation and could significantly decrease the 
number of candidates eventually detected. The red straight lines 
show the position of the two surveys if they were conducted over 
more pointings with the same total exposure time. Each dot rep- 
resent an additional pointing. 

we show the combination of area and limiting magnitude nec- 
essary to detect a given number of galaxies at 6.4 < z < 
7.2. We also show the position of our present data (labelled 
as Hawk-1 GOODS), of the planned Ultra- VISTA observations 
of the COSMOS fi eld, of the WFC3 pointing over the UDF 
jOesch et al.l l2009a) and of a ~20hr pointing with the 1 15 W fil- 
ter of JWST-NlRCam. For the Ultra- VISTA COSMOS survey, 
we consider the case of the deepest non-contiguous observations 
for a total of 0.73 deg^. We remark that these predictions do not 
include the fraction of galaxies lost because of different effects. 
In our case, the major sources of losses are the incompleteness 
in the photometric detection and the scatter (intrinsic and obser- 
vational) in the Z-Y colour that we use as threshold, leading to 
a loss of about 30% of the candidates. For this reason, the num- 
ber of candidates that we detect is about half of the expected 
numbers in Figure [15] Clearly, surveys with similar incomplete- 
ness levels should expect a comparable reduction of the observed 
numbers. The regression of the bright side of the LF, combined 
with its exponential slope, results in a flattening of the expected 
cumulative numbers at bright fluxes even in wide areas. In prac- 
tice, it is very difficult to detect galaxies brighter than 25 mag 
even over areas of about Ideg^, such that the detection of galax- 
ies brighter than 24.5 in Ultra- VISTA would be a formal viola- 
tion of our LF. 

The two red straight lines in Figure [15] show the expected 
position in the area-limiting magnitude plane of the two ground- 
based surveys (Hawk-1 GOODS and Ultra- Vista COSMOS) if 
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Fig. 16. Expected cumulative number counts of high redshift 
galaxies as a function of the F\6QW magnitude, for different 
redshift ranges and different choices of the evolutionary LF. The 
adopted colour cuts and the corresponding redshift ranges are 
shown in the legend. In both panels, the green dashed line is the 
prediction from a constant (from z = 6 to z = 10) UV LF. The 
black solid line is our best fit evolving LF, while the shaded area 
shows the predicted number counts for all the LF parameters ac- 
ceptable at the 1 cr level in our survey. The blu e continuous line 
is the evolving LF from lBouwens et al] (l2008l) . 



they were conducted over more pointings with the same total 
exposure time (32hrs and 150hrs, respectively). Increasing the 
number of pointings the total number of candidates decreases, 
with a loss of statistical robustness. For a relatively small number 
of pointings, this is likely counterbalanced by the decrease in the 
scatter due to cosmic variance. In essence, for a given instrument 
and total investment of time, the highest return in a statistical 
sense is given by a modest number of independent pointings. 

Only WFC3 on HST has the required IR sensitivity to de- 
tect galaxies at even higher redshifts, ahead of the advent of 
JWST. For this reason, several surveys are starting or are be- 
ing planned using the WFC3 IR filters. In Figure [T6]we provide 
the predicted cumulative number counts for galaxies in differ- 
ent redshift ranges, selected using either f 105M - F\25W > 
1 (that corresponds to the redshift range z - 7.8 - 9) or 
Film - F160W > 1 (z > 8.8). In both panels, the magni- 
tudes in the F\25W and F\6QW filters are computed from Ligoo 
using the same CB07 described above. We provide the number 
counts expected in the case of a non-evolving LF from z - 6 
(clearly an upper limit), the predictions from the evolving LF of 
Bouwens et al. (2008) (which was drawn from z = 4 to z = 10) 
and the extension to z = 10 of our LF. For the extrapolation of 
our LF we also show the whole region predicted extending from 
z = 7 to z = 10 all the LFs acceptable at the Icr level found from 
our analysis at z = 7. 

Apart from the unrealistic case of a constant LF beyond 
z = 6, the number density of LBGs at z > 8 is expected to 
be very low, and requires a combination of large areas and ex- 
treme depth to collect a sizeable sample. Taking into account the 
WFC3 size (about 4.6 arcmin^), it is necessary to reach contin- 
uum magnitudes at least mf\25w - w^Fieow - 28 to detect at least 
one dropout in a single WFC3 pointing. 

Figure [16] also shows that the evolution of the LF result- 
ing from this work is faster than the iBouwens et al.l (12008 ) one. 
Assuming that this evolution continues at the same rate at higher 
redshift, the number of expected dropouts would be smaller 
at higher redshifts. Clearly, deep and wide observations with 
WFC3 will be able to disentangle the various scenarios for the 
evolution of the UV LF. 
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8. Summary and Conclusions 

We present in this work the results of a F-band survey of the 
GOODS-South field, aimed at detecting galaxies at > 6.5 and 
measuring their number density. 

With this purpose we have made use of the deep 7-band ob- 
servations of the GOODS-South field obtained with Hawk-1, the 
new near-IR camera installed at the VLT. We have matched and 
combined these data with the publicly available images in the 
BVIZ (ACS), U and R (VIMOS) and JHK (ISAAC) bands. The 
final area covered by these observations is of about 90 arcmin^ 
at a magnitude limit Y ^ 26.7. 

We analyze this sample to select high-redshift (z > 6) candi- 
date galaxies following a Lyman-break colour criterion adapted 
to our filter set. Galaxies are selected in the Y band, and identi- 
fied by their large colour break Z-Y > 1 . Particular care has been 
taken in removing interlopers of various origins, including lower 
redshift galaxies with large Balmer breaks, variable sources and 
Galactic brown dwarfs, although some residual contamination 
from the latter cannot be excluded. An additional class of inter- 
lopers, exhibiting large Z-Y colours as well as significant emis- 
sion in the blue bands has been removed by requiring a strin- 
gent non-detection in the UBVRI bands. We argue that some of 
these contaminants might be very faint emission-line galaxies or 
AGNs at intermediate redshift. The accuracy and statistical im- 
pact of these criteria have been evaluated with extensive Monte 
Carlo imaging simulations. 

From the same simulations we estimate that our redshift se- 
lection function is mostly efficient in the interval 6.4 < z < 7.2. 

We eventually isolated 7 highly reliable z-drop candidates af- 
ter removing from the colour selected sample one known galac- 
tic cool dwarf star and one source undetected in deep available 
NICMOS F160W images. 

To estimate the constraints that our observations set on the 
evolution of the UV Luminosity Function at z > 6.5 we ran de- 
tailed and realistic imaging simulations of galaxy populations 
following different UV Schechter functions with linearly evolv- 
ing parameters log((p) and M, and a fixed value a - -1.71. Our 
simulations account for all the uncertainties involved in the ob- 
servations: detection completeness, photometric scatter and ran- 
dom fluctuations in the S/N measure due to overlapping unre- 
solved sources or other effects. We compare the resulting distri- 
butions of simulated magnitudes and colours with the observed 
ones following a maximum likelihood approach to constrain the 
parameters of the evolving UV LF. 

We find strong evidence for an evolution of the luminosity 
function above z=6: our analysis rules out at a ~ 99% confidence 
level that the LF remains constant in both <f> and M» above z - 6, 
even considering the effect of cosmic variance. 

From our maximum likelihood analysis we estimate - 
035t° \^^10-^Mpc-^ andM*(1500A) = -20.24 + 0.45 at flie me- 
dian redshift of our sample (z = 6.8). With respect to the values 
found at z = 6 our best fit model indicates an evolution both 
in the normalisation factor cp and in the characteristic magni- 
tude Mt. Our resu lts are consistent with the recent analysis by 
lOuchi etaTI (l2009l) . considering the large statistical errors and 
the high degeneracy between the two parameters. 

Fortunately, the uncertainty and the degeneracy in the M, 
and (f) best-fit values are not reflected in a comparable uncertainty 
in the number density of bright galaxies, i.e. in the integral of 
the LF up to M,. We derive a UV luminosity density puv - 
LS'll glO^^erg s"' //z"' Mpc^^ (M < -19) and a star formation 
rate density 5 = 2.8+f ^10"^M,„„ yr^' Mpc^l These values 



are definitely lower than the corresponding ones at z ~ 6 by a 
factor ~ 3.5. 

Although we did our best to carefully evaluate the systematic 
effects, we caution that our results depend on an extensive set of 
simulations to address the competing systematic effects which 
may increase or decrease the observed number of candidates. 
We note that the two most obvious - the possible residual con- 
tamination due to brown dwarfs and the possibility of a spurious 
detection for two of our candidates - will increase the amount 
of evolution from z = 6 to z = 6.8. Clearly, only spectroscopic 
surveys and/or deeper imaging in the IR will definitely settle the 
issue. 

Such a strong evolution in the UV LF has important conse- 
quences for reionization scenarios, as well as for the planning of 
future surveys aimed at detecting very high redshift galaxies and 
reionization sources. 

Determining whether the UV emission of normal galaxies is 
capable of reionizing the Universe at z> 6 requires the knowl- 
edge of the value of many parameters that are still unconstrained 
at these redshifts: the escape fraction of ionising photons, the 
Hll clumping factor, the exact spectrum of star forming galax- 
ies, their dust content as well as the shape of the stellar ini- 
tial mass fun ction and the metallicity of stellar populations (e.g . 
Ma dau et aL£ l999; Barkana & Loeb 2 ()0Tt IStiavelh et al.ll2004l) . 
Following B olton & Ha ehnelt (2007) and their conclusions on 
the early estimates from Bouwens et al. (2005) we note that the 
strong decrease we observe in the UV emission coming from 
relatively bright sources implies that this population alone is not 
capable of reionizing the Universe beyond redshift 6. These ob- 
servations can be reconciled with a completion of reionization 
before z = 6 only under the hypothesis of an evolution of the 
physical parameters quoted above, like an increase of the es- 
cape fraction, an harder UV sp ectrum, a lower dumpiness fac- 
tor or lower metallicities (e.g. iHenrv et"a]l 120091 : lOesch et"an 
2009b). Another possibility is that a relevant contribute to the 
UV emis sion comes f rom galaxies at the faint end of the LF 
(Bouwen s et al.l l2b07h or from more exo tic sources (see e.g. 
Yenkatesan et al.ll2003l:lMadau et al.ll2004l) . 

Thus it will be possible to fully assess the role of LBGs in the 
reionization of the neutral IGM only putting tighter constraints 
both on the bright and on the faint end of the luminosity function, 
at z~ 7 and beyond. 

Our results have also implications on the efficiency of these 
future surveys. We show that to collect sizeable samples of 
z > 7 bright galaxies it is necessary to reach very faint limits 
('WF160W S 28) over relatively large areas. On the other hand, 
the extrapolation of our results to higher redshifts indicates that 
dedicated surveys with HST-WFC3 will be able to disentangle 
the various scenarios for the evolution of the UV LF parameters 
above z ~ 6. 
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